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ABSTRACT 


The  basic  Michelson  interferometer,  as  used  for  Fourier  transform  spectroscopy,  is  analyzed  in  this 
report.  The  principles  of  operation  are  explained,  and  its  inherent  limitations  are  shown.  An  original 
analysis  of  apodization,  for  the  case  of  an  off-axis  detector  field-of-view,  is  included  in  an  appendix. 
Practical  limitations  of  real  instruments  are  also  analyzed.  These  include  misalignment,  wavefront  error, 
and  noise  sources.  In  another  appendix,  the  use  of  retroreflectors  in  the  interferometer  is  analyzed,  and 
the  effect  of  lateral  displacement  of  a  retroreflector  is  compared  to  that  of  mirror  tilt  in  a  plane-mirror 
interferometer.  Special  attention  is  given  to  sources  of  noise  in  the  Fourier  transform  spectroradiometer, 
and  expressions  are  derived  for  noise-equivalent  radiance  from  these  sources.  The  performance  of  a 
Fourier  transform  spectrometer  is  compared  with  that  of  a  filter-wheel  spectrometer  having  the  same 
optical  aperture  and  spectral  resolution.  The  mathematical  treatments  in  this  report  are  sufficiently  de¬ 
tailed  that  they  should  be  easy  to  follow,  once  the  optical  principles  are  grasped.  Some  numerical 
calculations  are  included,  based  on  the  specifications  of  the  proposed  GOES  High-Resolution  Interferom¬ 
eter  Sounder  (GHIS). 
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1.  INTRODUCTION 


Weather  satellites  use  infrared  spectrometers  to  measure  the  temperature  and  humidity  at  various 
levels  in  the  atmosphere,  in  a  process  known  as  atmospheric  sounding.  The  sounder  instrument  of  the 
newest  Geostationary  Operational  Environmental  Satellite,  GOES-8,  uses  a  filter  wheel  for  this  purpose. 
Lincoln  Laboratory  has  undertaken  a  study  of  the  feasibility  of  replacing  the  filter  wheel  on  future  GOES 
spacecraft  with  a  Fourier  Transform  Spectrometer  (FTS)  [1],[2].  This  report  presents  the  theory  of  Fourier 
Transform  Interferometry  (FTI),  as  it  relates  to  the  proposed  GOES  High-Resolution  Interferometer 
Sounder  (GHIS).  After  a  brief  introduction,  the  basic  theory  is  followed  with  a  discussion  of  the  perfor¬ 
mance  limitations  of  the  technique.  Finally,  the  performance  of  the  FTS  is  compared  with  that  of  a  filter 
wheel  system,  using  the  same  detectors  and  input  optics. 

Optical  spectra  were  first  measured  with  the  use  of  an  interferometer  by  Michelson  in  the  nine¬ 
teenth  century.  The  utility  of  the  technique  was  limited  at  the  time,  first  by  the  lack  of  sufficiently 
quantitative  detectors,  and  second  by  the  difficulty  of  performing  the  necessary  Fourier  transforms.  With 
the  advent  of  electro-optical  detectors,  and  especially  digital  computers,  those  barriers  were  removed  by 
1960.  Fourier  Transform  Spectroscopy  has  since  become  a  standard  tool  in  the  analytical  laboratory.  The 
Cooley-Tukey  Fast  Fourier  Transform  (FFT)  algorithm  (1965),  and  the  exponential  improvement  in  the 
cost/performance  ratio  of  computer  systems,  have  accelerated  the  trend. 

Many  Fourier  Transform  Spectrometer  (FTS)  instmments  have  been  flown  successfully  in  space, 
beginning  in  1962.  For  a  full  review  of  those  instruments,  see  the  report  by  Mooney,  et  al.  [3].  For  a 
history  of  Fourier  Transform  Spectroscopy,  and  a  detailed  treatment  of  its  principles,  see  Introductory 
Fourier  Transform  Spectroscopy  by  R.  J.  Bell  [4].  The  basic  theory,  and  instrumental  limitations,  are  also 
discussed  by  Brault  [5]. 
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2.  PRINCIPLE  OF  OPERATION 


The  simplest  type  of  FTS  is  based  on  a  Michelson  interferometer,  as  sketched  in  Figure  1.  The 
incoming  beam  of  collimated  radiation  is  divided  in  two  by  the  beamsplitter.  The  beam  in  each  leg  of 
the  interferometer  is  reflected  back  toward  the  beamsplitter  by  a  mirror.  At  the  beamsplitter,  each  beam 
is  again  divided,  one  part  going  to  the  detector,  the  other  returning  to  the  input  source.  In  each  direction, 
the  beams  from  the  two  legs  recombine,  or  interfere.  If  the  interferometer  mirrors  are  correctly  aligned, 
so  that  the  wavefronts  of  the  interfering  beams  are  parallel,  then  the  intensity  of  the  light  focused  by  the 
lens  on  the  detector  is  modulated  according  to  the  optical  phase  difference  of  the  beams  from  the  two 
interferometer  legs.  The  optical  phase  difference,  <(),  is  given  by  IkAxIX,  where  A  is  the  wavelength  of 
the  light,  and  Ax  is  the  optical  path  difference,  or  OPD.  If  yj  and  yj  are  the  optical  distances  of  the  two 
interferometer  mirrors  from  the  beamsplitter,  the  OPD  Ax  equals  2(y2  -  y^). 


Figure  1.  Optical  diagram  of  a  basic  Michelson  interferometer. 


In  operation,  one  or  both  of  the  interferometer  mirrors  are  moved  (while  maintaining  alignment) 
so  Ax  is  varied.  The  intensity  recorded  at  the  detector,  for  a  monochromatic  source,  is  then  given  by 


/^(v,  a:)  =  ^  /o[l  +  cos{27tvAx)\ . 


(1) 


3 


where 


V  =  MX  ,  the  wavenumber  (usually  measured  in  cm"*), 

Iq  =  the  intensity  at  zero  path  difference  (ZPD),  i.e.,  Ac  =  0. 

At  the  same  time,  the  beams  leaving  by  the  input  side  of  the  interferometer  interfere,  with  an 
intensity  given  by 


/,(v,jc)  =  i/o[l-cos(2m^Ac)]  . 


(2) 


Equation  (2)  is  given  here  merely  to  show  that  energy  is  indeed  conserved;  the  total  power  leaving  the 
interferometer  is  constant,  even  though  that  at  the  detector  is  modulated.  More  complex  instruments 
capture  both  outgoing  beams,  and  subtract  one  interference  signal  from  the  other.  However,  such  “bal¬ 
anced”  designs  are  not  considered  in  this  report,  since  they  require  a  second  detector  for  each  signal. 

In  the  foregoing,  it  was  assumed  that  both  the  reflectance,  /?,  and  transmission,  T,  of  the  beamsplitter 
were  exactly  one  half.  More  generally,  we  find 


2  tJ5[l  +  cos(27rvx)]  , 


(3) 


where 


Tj^  =  ART  =  the  beamsplitter  ejficiency.  (N.  b.:  0  <  7)^  <  1 .  77^  =  1  if  /?  =  T  =  0.5), 

=  the  transmission  of  the  remainder  of  the  optical  system,  excluding  the  beamsplitter. 

In  the  present  application,  the  input  is  a  spectral  radiance,  composed  of  many  wavenumbers, 
V,  the  total  signal  from  the  detector,  5(x),  becomes 

^00 

j(x)  =  2j^  5„[l-l-cos(2m'x)]dv  ,  (4) 

where  we  have  defined  a  spectral  signal,  as 

Sv  =  ^  ^  (5) 

and 

Ag  =  the  area  of  the  optical  aperture, 

£2  =  the  solid  angle  of  accepted  radiation,  and 

=  the  responsivity  of  the  detector  (in,  e.g.,  VoltsAVatt). 

[N.  b.;  The  product  often  called  the  throughput,  or  etendue,  of  the  system,  is  an  optical 
invariant.  It  is  unchanged,  from  the  input  aperture  to  the  detectors,  as  long  as  the  optics  are  properly 
designed  to  avoid  vignetting.] 
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The  detector  signal,  equation  (4),  is  composed  of  an  unmodulated,  or  “DC”  part,  and  a  modulated, 
or  “AC”  part,  of  average  value  zero.  The  value  of  the  unmodulated  component  is  just  equal  to  one  half 
the  signal  at  ZPD.  The  modulated  component  is  called  the  interferogram.  If  we  let  the  spectral  signal  be 
defined  symmetrically  for  both  positive  and  negative  wavenumbers  (i.e.,  S^,(-v)  =  S^,(v)),  then  the  inter¬ 
ferogram  can  be  expressed*  as 

m{x)  =  [  5y[cos(27ZVx)-t-isin(2m'x:)^v  ,or  (6) 

m{x)=r  .  (7) 

J— oo 

We  see  from  equation  (7)  that  m{x)  and  5^,  form  a  Fourier  transform  pair.  Hence,  the  spectral  signal 
could  be  recovered  from  the  interferogram  as 

5^  -  f”  m{x)e~^^^^dx  .  (8) 

J— oo 

Full  knowledge  of  the  other  factors  in  equation  (5),  many  of  which  depend  on  wavenumber,  then 
allows  us  to  determine  the  spectral  radiance,  The  very  important  subject  of  radiometric  calibration, 
however,  is  beyond  the  scope  of  the  present  report. 


*The  imaginary  term  in  the  integral  is  antisymmetric;  hence  it  vanishes. 
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3.  INHERENT  LIMITATIONS 


The  preceding  section  presented  the  basic  theory  of  the  interferometer.  Next,  we  describe  some 
inherent  limitations  of  the  Michelson  interferometer  as  a  spcctro.scopic  instrument. 

3.1  RESOLUTION 


The  interferogram.  mM,  is  what  is  directly  measured  by  recording  the  detector  output  as  a  function 
of  OPD.  However,  it  obviously  can  not  be  determined  for  x  from  to  oo!  In.stead,  there  is  some 
maximum  OPD,  L.  In  this  case,  .v  is  contained  in  [-L,  L],  since  we  use  a  two-sided  interferogram.  We 
obtain  an  approximation  to  the  spectral  signal  as 


5' =  I  ^m{x)e  . 

This  is  equivalent  to 


(9) 


where  FI  is  the  symmetric,  unit  rectangular  window  function.  Since  5,,  is  the  transform  of  a  product  of 
two  functions,  then  it  is  equal  to  the  convolution  of  their  (inverse)  transforms.  We  have 

5'.  =  2iL.sinc(2Z,v)*5,,  (the  symbol  *  stands  for  convolution).  (11) 

since  the  Fourier  transform  of  n(.v/2Z,)  is  given  by 


n 


2L 


=  2Lsinc(2Lv) : 


.sin(2niv) 


(12) 


From  equations  (II)  and  (12),  wc  find  that  the  effect  of  the  finite  interferometer  mirror  stroke  is 
to  limit  the  spectral  resolution  of  the  computed  spectrum.  The  true  spectrum  is  convolved  with  an 
instrument  resolution  function,  sinc(2Z,v),  which  has  its  first  zero  at 


By  convention,  8v  is  called  the  unapodized  spectral  resolution  of  the  interferometer. 


(13) 


3.2  APODIZATION 

So  far,  the  incoming  radiation  has  been  treated  as  having  a  negligible  angular  divergence  within 
the  interferometer  (i.e.,  £2  ~  0).  All  rays  were  assumed  to  be  parallel  to  the  interferometer  optical  axis, 
defined  as  the  normal  to  the  mirrors.  For  rays  at  an  angle  a  to  the  axis,  the  interferometer  operates  in 
the  same  way,  except  that  the  OPD  is  no  longer  x,  but  rather  .rcosa.  Returning  to  equation  (4),  the  detector 
signal  is  now  given  by 
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(14a) 


s{x)  =  2J  5',,  I cos(27rv'.r  cox  a)]sin  adad$  ijv 

’  e  I 


or 

•'■(•'■)  =  2 ‘^''[’■^^(“(>’^0)^'  (14b) 

where 


0  =  the  azimuthal  angle  about  the  interferometer  optical  axis,  and 


0  =  the  region  of  integration  representing  the  acceptance  ficld-of-view,  of  solid  angle  Q. 
and  we  have  defined  the  geometrical,  or  monochromatic  signal,  factor  in  equations  (14)  as 

G(ao,po)==—  f[cos(2Tuvxcosa)sinadou;l<t)  . 

Q  jJ 
0 


(15) 


The  integrand  in  equations  (14)  repre,sents  an  image  in  the  field  of  view  consisting  of  concentric 
circular  fringes.  They  have  an  appearance  similar  to  “Newton’s  rings,”  seen  when  a  spherical  optical 
surface  is  in  clo.se  proximity  to  a  planar  one.  The  rings  move  in  or  out  as  x  varies.  Figure  2  shows  a  series 
of  snapshots  of  the  fringes,  as  i’,r  ranges  from  1000  to  100  i  waves. 

Assuming  that  the  acceptance  field-of-view  is  circular,  as  usual,  but  that  it  may  be  displaced  from 
the  interferometer  optical  axis,  by  an  angle  a^,  we  can  simplify  the  limits  of  integration  by  changing  the 
angular  variables  to  p,  P,  as  defined  in  Figure  3. 

Equation  (15)  may  be  written  as 


G(ao,po)  =  — I  **  cos{2ffVx  cos  a) sin  pc(p<f/^  , 

and  the  solid  angle  of  acceptance  (included  in  the  definition  of  S^,)  is  given  by 

[f"  s'f  Pdpcip  =  2;r(l  -  cospo)  =  ;rp5  . 

In  case  a,j  =  0,  we  have  a  =p,  and  equation  (16)  becomes 
G((),  Po  )  =  ^  cos(2m'x  cos  p)  sin  pdp  . 

Making  the  substitution  y  =  cosp,  this  becomes 

G(0,  Po )  =  —  f  cos(2rtv>x:v)£/v  . 

Q  j  cos  Pi, 


(16) 


(17) 


(18) 


(19) 


from  which 


Figure  2.  Fringe  pattern  evolution  in  the  system  focal  plane,  as  the  on-axis  fringe  phase  varies  through  one  cycle. 
The  FOV  shown  is  30x30  mrad.  The  wavenumber  is  2700  cm'^,  and  the  interferometer  delay  starts  at  1000  waves 
on  the  axis. 


Using  equation  (17),  and  the  identity 


sin A-smB  =  2 sin  -  cos  - 

2  )  K  2  ) 


(21) 


we  arrive  at 


r 

i  -  cos  iTtvx 

1 — 

V  2  )  1, 

4n]) 

(22) 
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Figure  3.  Geometry  of  detector  acceptance  field-of-view  in  the  system  focal  plane.  Each  vector  represents  an  angle 
in  interferometer  space. 


If  we  use  the  definition  of  the  sine  function. 


sin  7DC 

sine  X  = - 

TlX 


equation  (22)  becomes 


(23) 


.  (vxQ^ 

\ 

sine  - 

cos 

lltvx 

\2n  ) 

\ 

4;rX 

Thus,  for  the  on- axis  field  stop,  equation  (14)  becomes 


(24) 


Comparison  of  equation  (25)  with  equation  (4)  reveals  the  consequences  of  having  a  non-zero  field-of- 
view,  £2.  One  effect  is  to  shift  the  apparent  wavenumber  (or  conversely,  the  apparent  OPD),  by  a  factor 
close  to  unity.  The  other  effect  is  to  modulate  the  fringe  amplitude  (or  contrast)  with  a  sine  function.  The 
argument  of  the  sine  function  contains  both  the  wavenumber  and  the  OPD.  The  reduction  of  the  fringe 
amplitude  as  bcl  increases  is  often  referred  to  as  self-apodization. 

The  exact  solution  of  the  integration  problem  in  the  more  general  case,  that  ttg  0,  may  be  found 
in  Appendix  A,  in  terms  of  Lommel  functions.  Again,  as  a  result  of  the  non-zero  aperture  diameter,  the 
interferogram  is  modulated  by  an  overall  factor,  M,  which  reduces  the  fringe  contrast  as  the  OPD 
increases: 


'  .  (vxQ\ 

r, 

1  +  sine  -  cos  2jtvx 

1 — 

yiK )  y 

4;r. 

JJ 

dv 


(25) 


G(ao.Po)“M«o.Po)cos 


iTtvx  cos  ag 


\  sin^Po^ 
4 


(A-34) 


10 


M  is  well  approximated,  for  at  least  1,000  waves  of  OPD,  by 


M«0>Po)  = 


,  {l7Wx)^Q(  .  2  ^2 

l--i — 1  sin^  oco  +  — 


%7l 


VlK 


(A-35b) 


In  addition  to  the  fringe  visibility  reduction,  the  apparent  wavenumber  is  altered,  by  a  factor  which 
depends  on  both  the  aperture  size  and  its  offset  angle: 


v'  =  vcosag 


\1l 


(A-36) 


By  expanding  the  sine  function  as  a  power  series,  the  reader  may  easily  verify  that  equation  (A-34) 
reduces  to  the  same  solution  as  in  equation  (24),  when  goes  to  0. 


3.3  THROUGHPUT 


The  preceding  section  revealed  that  there  is  a  trade-off  to  be  made  between  optical  throughput, 
A.qQ.,  and  spectral  resolution,  8v.  This  comes  about  because  of  the  solid  angle  factor  in  the  apodization 
function  in  either  equation  (24)  or  equation  (34).  The  modulation  in  the  interferogram  approaches  zero 
as  the  product  vxQ  approaches  1%.  The  maximum  absolute  value  of  x  is  L,  which  gives  the  spectral 
resolution  (see  equation  (13)).  Writing  the  sine  function  in  terms  of  spectral  resolving  power, 


we  have 


Psr^ 


V 

5v 


sine 


2n  j 


=  sine 


Psr^) 

471  J 


(26) 


Assuming  that  the  application  calls  for  a  certain  resolving  power  at  a  certain  wavelength,  one 
normally  wishes  to  use  the  maximum  solid  angle  to  get  the  best  signal-to-noise  ratio.  Since  the  interfero¬ 
gram  signal  (see  equations  (4)  through  (7))  also  includes  the  solid  angle  factor,  we  wish  to  maximize  the 
product 


^3slnc  — —  = - sin 


\  Alt  )  P,,  V  4 

Staying  within  the  central  lobe  of  the  sine  function,  the  maximum  signal  is  obtained  for 

lit 


(27) 


42 


MAX 


(28) 


Thus,  there  is  an  inverse  relationship  between  the  optical  throughput,  which  is  proportional  to  solid  angle, 
and  resolving  power. 
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4.  PRACTICAL  LIMITATIONS 


The  preceding  discussion  concerned  the  theoretical  capabilities  of  the  Michelson  interferometer. 
Next,  we  turn  to  the  performance  limitations  of  practical  hardware,  caused  by  the  effects  of  optical 
misalignment,  fabrication  tolerances,  and  noise. 

4.1  MISALIGNMENT 

So  far,  we  have  implicitly  assumed  that  the  interferometer  is  perfectly  aligned,  that  is  to  say,  each 
mirror  is  exactly  parallel  to  the  image  of  the  other  mirror,  formed  by  the  beamsplitter.  Suppose  now  that 
one  mirror  is  misaligned  by  a  small  angle,  £.  The  optical  path  difference  is  now  a  function  of  the  position 
of  the  rays  within  the  aperture.  Let  the  polar  coordinates  (r,  (p)  denote  the  position  of  the  rays  in  the  plane 
of  the  Lyot  stop,*  or  aperture  of  the  system.  For  these  rays,  the  optical  path  difference,  x,  in  our  original 
expressions  is  replaced  by 

x'  =  x  +  2ersin<j)  .  (29) 

The  total  signal,  given  by  equation  (4),  now  becomes 


s(j£:)  =  2j^  5^,[l  +  //£(v,x)]^iv 

where 


(30) 


cos[2m'(x  +  2er  sin  ^)]  rdrd<j>  ,  and 
R  =  the  radius  of  the  Lyot  stop. 


To  evaluate  we  first  expand  the  cosine; 

kk 

where 


(31) 


(32) 


Fj  =  cos(2m^x) J  cos{47Cversin^)rdrd(l>  ,  and 

F2  =  sin(2;rvx)J  sin(4;rv£rsin^)r<in/^  =  0  . 

Making  use  of  the  fact  that  [7] 

J  cos[xsm.<j))d^  =  7cJQ{x)  ,  and  j x'^J„_i{x)dx  =  x'‘J„{x) 


(33) 

(34) 


*The  Lyot  stop  is  a  defining  aperture  employed  to  minimize  stray  light  reaching  the  detector.  It  was 
introduced  in  a  solar  coronagraph  developed  by  B.  Lyot  in  1931.  (See  R.  Klingslake,  Optical  System 
Design  (1983),  pp.  228-229  [6].)  In  an  FTS,  it  is  frequently  located  at  the  output  of  the  interferometer. 
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leads  to 


H,{v,x)  = 


2 

—^cos 

FT 


{27CVx)j  jQ{47CV£r)rdr 


from  which  follows 


(35) 


where 


Hg{v,x)  =  —  Ji{k)cos{2m’x) 
k 


(36) 


k  =  AtcveR  ,  (37) 

We  thus  find  that  the  interferogram  is  multiplied  by  a  function  of  exactly  the  form  found  for 
diffraction  from  a  circular  aperture.  For  k  =  0,  its  value  is  exactly  1.  The  first  zero  occurs  for  k  =  3.83. 
For  ^  «  1,  we  can  use  a  power  series  expansion,  retaining  only  the  first  two  terms: 


1  (k 


-jAk)  =  \--  -  +  — 

I-*'-'  I-*! 


2(2 
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\-2{7tx’eRy 


The  condition  k  «  I  amounts  to 


(38) 


£« 


A 

47tR 


(39) 


which  implies  that  the  linear  displacement  of  the  mirror  must  vary  by  «  X/2n:  across  its  diameter. 
Otherwise,  the  fringe  contrast  will  be  significantly  reduced. 

In  view  of  the  sensitivity  of  the  Michelson  interferometer  to  tilt  misalignment  of  the  flat  mirrors, 
one  might  use  comer-cube  retroreflector  mirrors.  However,  retroreflectors  also  have  a  sensitivity.  In  their 
case,  it  is  to  translational  alignment,  as  explained  in  Appendix  C. 


4.2  WAVEFRONT  ERRORS 


Michelson  interferometers  are  quite  sensitive  to  anything  which  causes  a  difference  in  the  wavefronts 
of  the  two  beams  when  they  are  recombined.  This  has  just  been  shown  for  the  case  of  misalignment.  We 
now  consider  the  effects  of  wavefront  distortions,  or  optical  path  difference  (OPD),  caused  by  departures 
from  flatness  in  the  beamsplitter  and  interferometer  mirror  surfaces.  In  this  treatment,  the  nominal  OPD, 
X,  caused  by  displacement  of  the  moving  mirror,  will  be  referred  to  as  the  optical  delay.  The  OPD  from 
component  distortions  will  be  called  wavefront  error.  This  wavefront  error  is  actually  the  difference 
between  the  wavefronts  from  the  two  interferometer  arms,  after  recombination  at  the  beamsplitter,  and 
after  removing  the  overall  optical  delay. 

The  analysis  of  wavefront  errors  begins  almost  exactly  like  that  of  misalignment  (indeed,  misalign¬ 
ment  can  be  viewed  as  a  special  case  of  wavefront  error): 
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5(x)  =  2j  5,,[l  +  7/2(v,Jc)]rfv  , 


where 


/f^(v,;r)  =  -^2  Iq ‘^^2;rv(x  +  5(r,^))]r^/^•^/^  ,  (41) 

R  =  the  radius  of  the  Lyot  stop, 

jc  =  optical  delay,  between  best-fit,  flat  wavefronts  from  the  two  arms,  and 
8(r,<j))  =  wavefront  error  at  radius  r,  azimuth  in  the  Lyot  stop. 

Evaluation  of  the  integral  over  the  Lyot  stop  begins  as  in  the  case  of  tilt  misalignment,  starting  with 
an  expansion  of  the  cosine: 

=  >  (42) 

kR 

where 

Pj  =  cos(2;rvx) J  cos{27CvS[r,<j>))rdrd(j>  ,and  (43a) 

P2  =  sin(2m'x) J  sm{27uv5{r,<t))ydrd<j)  .  (43b) 

The  quantity  v5(r,0)  is  the  wavefront  error,  expressed  in  wavelengths,  at  wavenumber  v.  We  know 
that  this  should  be  much  less  than  one  wave,  in  order  to  obtain  an  effective  interferogram,  so  we  use 
power  series  expansions  of  the  cosine  and  sine: 

1  2 

cos(2m'5(r,^))  =  l-— (2m’5(r,0))  H —  ,and  (44a) 

1  3 

sin(27Cvd[r,(pyj  =  27tvS[r,(l))-—^27CvS{r,(l))^  H —  .  (44b) 

We  choose  to  define  5  such  that  its  mean  value  over  the  aperture  is  zero: 

I  d{r,(l))rdrd(t)  =  0  .  (45) 

Also,  we  make  the  approximation  that  the  third  and  higher  powers  of  v5  average  to  zero  over  the  Lyot 
stop.  Making  the  substitutions  in  (43a),  we  have 


where 


Pj  »  cos(2m'x:)  ;rP^  -  2;r^  J"  J"^  (v5(r,^)]  rdrd(j)  ,or 
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(47) 


Pi  ~  tcR^  cos(2m^j!:)|l  —  'liP'  , 

where 

the  mean-squared  value  of  the  wavefront  error. 

Thus  we  find 

Hs{v,x)^cos{27tvxj^\-2n^v^{5'^)]^  .  (48) 

Once  again,  the  interferogram  contrast  is  now  reduced  by  a  modulation  factor,  given  by 

Mg^\-27t'^v^l^5'^)  ,  (49) 

or  more  simply,  the  fringe  contrast  is  reduced  from  its  ideal  value  of  one  by 

-AMg^2n^v^(^S^'^  ,  (50) 

The  total  wavefront  error  over  the  aperture  is  related  to  errors  introduced  by  the  various  compo¬ 
nents  (see  Appendix  B)  as 

^(r,  0)  =  2[5^5-5;i^(;' -I- 5^7  +  253  cos  0]  ,  (B-4) 

where 

Sju^ir,  0)  =  the  surface  irregularity  of  the  mirror  in  the  splitter  arm, 

Sf^c^r,  0)  =  the  surface  irregularity  of  the  mirror  in  the  compensator  arm, 

Sjfj-,  0)  =  the  surface  irregularity  of  the  beam-splitting  surface, 

0)  =  the  transmitted  wavefront  error  of  the  beamsplitter, 

0)  =  the  transmitted  wavefront  error  of  the  compensator,  and 
0  =  the  angle  of  incidence  at  the  beamsplitter. 

Appendix  B  gives  further  expressions  which  will  be  useful  for  estimating  the  loss  of  fringe  contrast 
caused  by  specified  tolerances  on  the  surface  irregularities  and  transmitted  wavefront  errors  of  the  various 
optical  components.  Here,  it  is  worthwhile  to  make  two  observations:  First,  irregularities  of  the  splitting 
surface  (Sg)  cause  roughly  twice  as  much  error  as  those  on  either  mirror,  because  a  negative  deviation 
in  one  arm  of  the  interferometer  is  subtracted  from  a  positive  deviation  in  the  other.  Second,  identical 
irregularities  on  the  two  mirrors,  or  identical  errors  on  the  two  transmissive  parts,  will  cancel  out. 
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4.3  SPLITTER-COMPENSATOR  DIFFERENCES 


A  compensator  plate  was  alluded  to  in  the  preceding  subsection,  but  its  function  has  not  yet  been 
explained.  It  is  needed  to  balance  the  optical  path  length  introduced  by  the  plate  (or  substrate)  on  which 
the  beam-splitting  surface  is  formed.  If  the  index  of  refraction  of  the  splitter  substrate  were  constant,  no 
compensation  would  be  needed,  as  the  shift  of  one  mirror  would  have  the  same  effect.  However,  the 
phenomenon  of  dispersion  in  the  substrate  requires  the  addition  of  a  compensator  of  exactly  the  same 
composition  and  thickness  as  the  splitter  substrate. 

There  must  be  some  tolerance  on  the  thicknesses  of  the  beamsplitter  and  compensator  plates,  as 
actually  fabricated.  When  a  plate  of  thickness  t  is  introduced  into  a  beam  at  angle  of  incidence  0,  the 
optical  path  is  increased  by 


A  =  t[n  cos  6'  -  cos  0]  =  r 


i 


-  sin^  0  -  cos  0 


(51) 


The  difference  in  optical  paths,  between  the  compensator  and  splitter  plates,  when  their  thicknesses  differ 
by  At,  is 


S  =  At] 


■Jn^ 


sin"^  0  -COS0 


(52) 


The  change  in  this  OPD,  between  wavenumbers  Aj  and  where  the  indices  of  refraction  are  «j  and  n2, 
respectively,  is 


=  At  -Jn^  -sin^  0  -  -sin^0 


(53) 


This  may  be  written  as 


1  A  2 

jnf  -  sin^  0 

j 

ll-f  2  9  1 

nf  -  sin^  0 

where  we  define 


(54) 


A  2  2  2 

An  s  «2  —  «! 

Assuming  that  An  «  n,  the  bracketed  quantity  in  (54)  may  be  expanded  as  a  binomial  series  to  yield  the 
approximation 

_  AtAn^ 

The  beamsplitter  and  compensator  for  GHIS  are  to  be  made  of  KBr.  The  dispersion  equation  for 
the  index  of  refraction  of  KBr  at  22°C  is  (see  [8],  pp.  7-70*) 


*Typographical  errors  in  the  formula  as  printed  in  the  handbook  are  corrected  in  this  report. 
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2  ^  i  r>-4  i2  i/^-8'i4  7.676x10  1.56569x10 

n  =2.361323-3.11497x10  A  -5.861x10  A  + - t ^ - 

A^  A^- 0.0324 

where  the  wavelength  is  given  in  /tm.  In  the  wavelength  range  of  GHIS  (3.7  ~  15  jum),  we  find  that 
-3.11x10“^  (a? . 

With  6  =  30°  and  Aj  =  15  ^m,  we  have 


:  0.0245 


^  Ji2 
15fJm 


At 


The  effect  of  the  dispersive  offset  of  the  optical  paths  in  the  two  interferometer  arms  is  to  introduce 
an  asymmetry  in  the  interferogram.  The  zero  path  difference  position  is  truly  zero  only  at  one  wavenumber; 
at  other  wavenumbers,  it  is  shifted.  When  the  Fourier  transform  is  performed  to  obtain  the  raw  spectrum 
(see  equation  (9)),  an  asymmetry  will  cause  the  appearance  of  an  imaginary  component  in  S^.  The  raw 
spectrum  will  contain  a  complex  phase  factor,  of  phase  The  calibration  process,  carried  out  in 

either  the  interferogram  or  complex  Fourier  transform  domains,  should  correct  for  this  shift.  However, 
as  a  guideline  for  system  design,  it  is  suggested  that  this  phase  angle  be  held  to  less  than  Itt,  in  order 
to  prevent  ambiguities  in  recovery  of  the  phase. 


4.4  NOISE  SOURCES 


The  interferogram  measurements  from  a  real  instrument  will  be  affected  by  noise  from  many 
potential  sources.  First,  we  shall  examine  the  relationship  between  the  noise  in  the  real  interferogram  and 
the  resulting  apparent  noise  in  the  computed  spectrum.  Then,  the  important  sources  of  noise  will  be 
considered  in  turn. 


4.4.1  Conversion  from  Interferogram  to  Spectral  Domain 

In  practice,  the  interferogram  is  measured  as  a  series  of  discrete  samples,  taken  at  N  different 
points.  As  before,  we  assume  a  two-sided  interferogram,  so  that  -L  <  x  <  L,  and  the  spacing  between 
samples  is  thus 


(56) 


A  Discrete  Fourier  Transform  is  performed  on  the  interferogram,  fx),  which  gives  a  computed 
spectrum,  F(v),  also  having  N  points.  The  spectrum  F{y)  is  periodic,  with  period  2v^.  The  following 
relations  hold  [5], [9]: 


N  = 


dx 


2v, 


5v 


(57) 
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We  shall  make  use  of  the  Rayleigh  Power  Theorem,  whieh  in  discrete  form  is  stated  as 


=  ■  (58) 

m=0  j=0 

The  theorem  may  also  be  restated  as 

2l(|/P)  =  2v„(|F|2)  ,  (59) 

where  the  brackets  denote  the  mean  value  over  the  whole  period  of  the  function. 

In  particular,  suppose  that  in  place  of /(a:),  we  put  a  real  measurement  noise,  in  the  interferogram 
domain.  We  then  have  a  direct  relationship  for  the  resulting  noise  in  the  spectral  domain, 

cr^V^  =  <Tvi/2v„  .  (60) 

From  the  relations  (57)  above,  we  find 


(61) 


At  this  point,  we  make  use  of  the  fact  that  because  the  spectrum  of  incident  light  must  be  real,  its 
true  interferogram  is  symmetric.  There  is  no  such  constraint  on  most  sources  of  noise.  Generally,  the  noise 
power  in  the  interferogram  will  be  evenly  divided  between  symmetric  and  antisymmetric  parts,  which 
correspond  to  real  and  imaginary  noise  components,  respectively,  in  the  computed  spectrum.  The  imagi¬ 
nary  half  of  the  noise  power  can  be  disregarded,  so  that  the  rms  (root-mean-square)  noise  in  the  spectrum 
is  divided  by  V2; 


4.4.2  Noise-Equivalent  Spectral  Radiance 


(62) 


We  specify  the  spectroradiometric  performance  of  the  instrument  in  terms  of  a  noise-equivalent 
spectral  radiance,  denoted  by  NEN^.  This  is  just  a  radiance  which  would  produce  an  output  signal 
equivalent  to  the  output  noise,  at  each  apparent  wavenumber.  If  cr^  is  the  system  spectral  noise,  then  the 
noise-equivalent  spectral  radiance  is  given  by 


NEN,=^B,  . 

Substituting  equation  (5)  for  the  signal,  we  have 


NEN^ 


4gv 


(63) 


(64) 
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The  system  NEN^  is  composed  of  contributions  from  various  sources.  Some  of  them  are  multipli¬ 
cative,  i.e.,  proportional  to  the  input  radiance.  Most,  however,  are  additive  noise  effects  in  the  interfero- 
gram.  Multiplicative  noise  sources,  which  include  interferometer  tilt  jitter  noise  and  velocity  error  noise, 
are  evaluated  as  noise-to-signal  ratios.  The  total  of  additive  errors  in  the  interferogram  can  be  treated  as 
the  sum  of  the  noise  powers  from  the  various  individual  sources: 

^  a^l  ,  {/}  =  sources  of  noise. 
i 

Each  term  in  the  summation  corresponds  to  a  noise-equivalent  spectral  radiance  contributed  by  that 
source: 


i 

We  evaluate  each  noise  contribution  in  terms  of  NEN^  in  order  to  have  a  common  basis  for 
comparison  of  the  different  sources  with  each  other,  and  of  the  system  NEN^  with  the  required  perfor¬ 
mance  specification. 

Upon  expressing  the  spectral  noise  (equation  (64))  in  terms  of  the  real  noise  in  the  interferogram 
(equation  (62)),  we  find 


NEN^ 


[Y 

V  N 


4.4.3  Photon  and  Detector  Noise 


(67) 


A  very  important  class  of  additive  noise  appears  at  the  output  of  the  detector.  It  is  made  up  of 
photon  noise  and  detector  noise.  The  photon  noise  (also  called  quantum  noise)  generally  has  a  uniform, 
or  “white”  temporal  spectrum,  and  it  results  from  the  inevitable  statistical  fluctuations  in  the  number  of 
photons  converted  to  photoelectrons.  (We  assume  that  photon  detectors  are  used  here.)  The  detector  noise 
can  have  several  components.  In  semiconductor  infrared  detectors,  charge  carriers  are  generated  by 
thermal  excitation  across  the  band  gap.  In  photoconductive  detectors,  this  is  referred  to  as  “thermal 
generation-recombination”  noise,  or  “thermal  G-R”  noise.  Photovoltaic  detectors  have  a  steady  “dark 
current,”  which  causes  noise  through  its  statistical  fluctuations. 

Both  photon  noise  and  intrinsic  detector  noise  are  usually  lumped  together  and  treated  as  the  total 
“detector  noise,”  which  we  shall  designate  When  divided  by  the  detector  responsivity,  it  gives  the 
detector’s  noise-equivalent  power: 


^x,det 


(68) 


The  usual  figure-of-merit  for  specifying  detectors  is  the  specific  detectivity,  or  D*.  It  is  the  recip¬ 
rocal  of  the  NEP,  but  normalized  to  unit  detector  area  and  temporal  bandwidth: 
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(69) 


NEPdet 


where 

=  the  detector  area,  and 
Af  =  the  electronic  bandwidth. 

(Notice  that  the  conditions  of  the  photon  flux  incident  on  the  detector  must  be  specified,  along  with  D*. 
In  the  present  work,  we  calculate  D*  based  on  the  photon  irradianee  expected  in  use,  and  not  on  the  more 
typical  295  K  blackbody  radiance  over  the  full  hemisphere.) 

The  effective  noise  bandwidth,  Af,  in  this  sampled-data  system  is  equal  to  the  Nyquist  frequency, 
which  is  one-half  the  sampling  frequency: 


where 


Af  fNy  2  fsamp 


N_ 

2T  ’ 


(70) 


T  =  the  duration  of  the  interferometer  scan,  of  N  samples. 

Returning  to  the  expression  for  NEN^  in  terms  of  interferogram  noise  (equation  (67)),  and  using  the 
foregoing  definitions,  we  find  a  useful  expression  for  the  detector  NEN^: 


NEN, 


4L 


v,det 


A^QXoxrisD* 


(71) 


Using  the  definition  of  the  unapodized  spectral  resolution,  5v,  this  may  also  be  written 


NEN, 


V,det  ~ 


AoEiToxlsSvD* 


4.4.4  Electronics  Noise 


(72) 


The  electronics  following  the  detector  also  contributes  noise.  It  mainly  originates  in  the  input  stages 
of  the  preamplifier.  It  may  be  characterized  by  a  noise  power  spectral  density,  which  is  referred  to 
the  input  of  the  preamplifier,  i.e.,  the  output  terminal  of  the  detector.  It  contributes  an  rms  noise  to  the 
interferogram  given  by 

^x,amp  ~  ^amp  V Af  ~  ^amp  '  (7^) 

The  NEN^  component  contributed  by  electronics  noise  is  thus 
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(74) 


MEN., 


4Ln 


amp 


4.4.5  Aliased  Noise 


When  the  discrete  Fourier  transform  of  a  sampled  signal  is  performed,  only  frequency  components 
below  the  Nyquist  limit  are  faithfully  preserved.  Any  signal  component  beyond  the  Nyquist  frequency 
appears  in  the  spectmm  at  its  alias  frequency,  given  by 

f alias  =  fsamp  "  /  =  '^fNy  “  /  (assuming  that  fNy<f<  2fMy)-  (75) 

The  raw  signal  must  be  filtered  to  eliminate  components  above  the  Nyquist  frequency  before 
sampling.  Once  they  are  sampled,  such  aliased  components  are  inseparable  from  the  true  signal.  Such 
filtering  is  done  first  in  the  optical  domain,  by  passband  filters  and/or  dichroic  beamsplitters.  Next,  there 
is  the  spectral  response  characteristic  of  the  detector.  Finally,  the  electrical  output  of  the  detector  is 
filtered.  In  the  interferometer,  the  temporal  frequency  is  directly  related  to  the  wavenumber  by 

f  =  v—  =  vVoPD  ,  (76) 

where 

Vqpd  ~  optical  path  difference  velocity. 

Heretofore,  the  electronic  filter  was  assumed  to  have  no  effect  on  the  signals  being  measured.  In 
other  words,  its  magnitude,  Mif),  was  taken  as  1  exactly.  To  be  correct,  we  must  substitute  M(f)dl^ 
wherever  appears  in  the  equations  above. 

Aliasing  will  produce  a  noise  component  corresponding  to  the  combined  detector  and  electronic 
rms  noise  amplitude  times  the  electronic  filter  magnitude  at  the  image  frequency: 

alias  =  ^{fimage)^^EN^,del  +  ^^N^amp  ’OT  (77) 


where 


alias 


4LM(f image) 


xVs^{^^OPd)  ^ 


•  + 


*^amp 

V  y 


(78) 


fimage  ^^OPD  •  (79) 

4.4.6  Low-Frequency  (“1//”)  Noise 

Semiconductor  devices,  including  photoconductive  IR  detectors  in  particular,  generally  exhibit 
excess  low-frequency  noise.  It  is  often  called  “1//”  noise,  because  its  power  spectral  density  has  the  form 

ayj  =  Kf-‘^  ,  (80) 
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where  a  is  a  constant  on  the  order  of  unity. 

In  device  specifications,  the  magnitude  of  the  Mf  noise  is  usually  given  in  terms  of  the  “white” 
noise,  n^,  at  higher  frequencies,  by  specifying  the  “knee  frequency,”  at  which  the  two  components 
are  equal.  The  total  noise  power  takes  the  form 


2  2 
n  =  «o 


^  ^  A 

Jknee 


1  + 


/ 


(81) 


The  1//  component  of  noise  from  the  detector  and  amplifier  can  thus  be  included  in  the  D*  and 
NEN^  expressions.  Notice  that  the  frequency  dependence  of  these  noise  components  now  implies  a 
wavenumber  dependence  in  NEN^  (see  equation  (76)). 

We  saw  earlier  that  electronic  filtering  is  needed  to  attenuate  high-frequency  noise,  so  that  it  would 
not  appear  on  the  signal  through  aliasing.  A  high-pass  filtering  should  be  used  as  well,  to  attenuate  the 
l/f  noise,  at  frequencies  below  the  low-wavenumber  end  of  the  signal  spectmm. 


4.4.7  Quantization  Noise 


The  analog-to-digital  converter  (ADC)  introduces  an  rms  error  in  the  digitized  signal,  which  is 
given  by 


^x,ADC 


2  •  w(0) 


(82) 


where 

m(0)  =  the  maximum  anticipated  interferogram  signal,  at  the  ZPD  peak,  and 

^nom  ~  l-^  =  the  effective  number  of  bits,  accounting  for  ADC  nonlinearity  errors. 
Referring  to  equations  (5)  and  (7),  we  can  write 


m{0)  =  ^A„n\  To^ri,9l^M{vVoPD)Bydv  . 

A  J — oo 


The  rms  ADC  error  is  thus 


_  o 


A^nNd\^ 


*  x,ADC  l) 


2-2^^ 


(83) 


(84) 


where  we  have  defined 


{roxils%M{vVoPD)By)  = 


2 

N5v 


r„xns%^{vyoPD)Bvdv  . 


(85) 
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When  is  substituted  in  the  general  NEN^  expression  (equation  (67)),  we  find 


MEN, 


1 


v,ADC  ~  h 


i^ox^s^^^iy^OPE^^v)  I N 


'^ox'nsK^iv^OPD)  V  6 


(86) 


4.4.8  Mirror  Velocity  Errors 


Ideally,  the  interferometer  signal  is  sampled  at  uniform  increments  of  optical  path  difference.  Ax. 
To  this  end,  a  reference  interferometer  provides  convenient  position  markers  to  trigger  the  sampling.  The 
mirror  velocity  does  not  affect  the  measurement  directly.  However,  the  temporal  frequency  of  the  signals 
is  proportional  to  the  mirror  velocity.  Since  the  overall  transfer  function  of  the  detector  and  electronics 
is,  in  general,  not  constant  with  frequency,  a  change  of  velocity  will  affect  the  measured  amplitude  and 
phase  of  any  given  spectral  component. 

We  refer  to  the  modulation  in  the  magnitude  of  the  spectral  signals  as  mirror  velocity  error.  It  is 
characterized  by  a  Noise-to-Signal  Ratio  (NSR)  given  by 


NSR, 


v,VEL  -  ■ 


^^VEL 

1  dM 

V2 

M{f)  df 

/=vV, 


OPD 


where 

^VEL  ~  variation  in  the  OPD  velocity,  and 

M(f)  =  the  magnitude  of  the  overall  transfer  function  at  frequency  /. 


(87) 


4.4.9  Sampling  Errors 

It  is  to  be  expected  that  the  signals  which  trigger  the  sampling  will  contain  some  position  error, 
whatever  their  source.  An  error  in  the  position  of  the  sample  produces  an  error  in  the  measured  amplitude 
of  the  interferogram,  which,  for  any  given  wavenumber,  varies  as  the  cosine  of  the  OPD.  The  case  of 
random  (uncorrelated)  sampling  errors  has  been  analyzed  by  Bell  and  Sanderson  [10].  Expressed  as  a 
noise-equivalent  radiance,  the  effect  is 

NENsAMP  =  27td[vB{v)]^^^  ,  (88) 

where 

•  (89) 

(vl,  v2)  =  the  range  of  detectable  wavenumbers,  and 
8  =  the  rms  position  error  in  the  sampling. 
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The  position  error,  5,  here  includes  error  contributions  from  noise  in  the  reference  interferometer, 
and  from  electrical  phase  shifts  induced  by  mirror  velocity  variations.  The  reference  interferometer  pulses 
will  have  an  equivalent  position  error  given  by 


^x,REF 


^REF 

7,7t  •  SNRj^£^p 


(90) 


where 

^REF  ~  reference  wavelength,  and 

SNRj^^f  =  the  signal-to-noise  ratio  from  the  reference  detector  and  amplifier. 

Variations  in  the  OPD  velocity  will  shift  the  electrical  frequency  of  the  interferometer  signals,  as 
described  earlier.  This  affects  the  phase,  0,  of  the  electrical  transfer  function,  as  well  as  the  magnitude. 
The  group  delay  of  the  signal  is 

_ _ 1_  ^ 

Itc  df 

Both  the  signal  being  measured  and  the  reference  signal  have  group  delays,  and  it  is  the  difference  of 
these  delays,  times  the  random  OPD  velocity  variation,  which  produces  a  random  position  error.  When 
combined  with  the  reference  position  error,  we  have,  for  the  total  sampling  position  error, 

5  =  ^ (cTypi (Tffi  -  +  (jisEF  .  (92) 

In  the  foregoing,  the  velocity  variations  were  assumed  to  be  random,  i.e.,  uncorrelated  from  sample 
to  sample.  In  practice,  the  velocity  will  probably  vary  slowly,  on  the  time  scale  of  the  sampling.  A  narrow 
spectrum  of  velocity  variations  will  result  in  relatively  smaller  sampling  errors  [11].  Thus,  the  expression 
given  here  for  f^E^sAMP  (equation  (87))  may  be  regarded  as  an  upper  limit. 

4.4.10  Optical  Jitter-Induced  Noise 

In  Section  4.1,  it  was  shown  that  the  amplitude  of  the  interferogram  is  reduced  by  any  misalignment 
of  the  interferometer.  The  modulation  factor  is 

Mie)  =  hi{k)^l-2{KveRf  ,  (93) 

/C 

where  as  before 

e  =  tilt  error  of  interferometer  mirror,  and 
R  =  radius  of  Lyot  stop. 

The  tilt  error,  e,  will  be  composed  of  both  a  systematic,  or  repeatable  tilt,  Eq,  and  a  random  tilt,  a, 
also  called  optical  jitter: 


\f=vVc 


OPD 


(91) 
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£  =  eg  +  a  . 


(94) 


To  find  the  modulation  caused  by  the  jitter  alone,  we  perform  a  Taylor  series  expansion  about  the 
repeatable  tilt  position,  Eq-. 

AM  = -2{m>Rf‘^eQa  +  .  (95) 

The  variance  of  M  is  given  by 

ali^i^{AMf)-{AMf  .  (96) 

In  terms  of  the  variance  and  mean  of  the  tilt  components,  this  is 

ah  =  A{TcvRf  Aiel  )al  +  2a^  ,  (97) 

where 

{a)  =  0  al^(a^)  . 

The  modulation  caused  by  the  jitter  results  in  a  multiplicative  noise  in  the  spectrum,  given  by 

NSRjiyt  =  =  2{7tvRf  a  a  ^2^£o)  +  0’a  •  (98) 

The  square  root  of  two  reduction  in  the  noise  comes  from  the  fact  that  only  the  real  noise  power,  and 
not  the  imaginary  component,  contributes  to  noise  in  the  reconstructed  spectrum,  for  the  case  of  two-sided 
interferograms. 

4.4.11  Total  Noise-Equivalent  Spectral  Radiance 

To  summarize  the  foregoing  results,  the  total  noise-equivalent  spectral  radiance  is  given  by 
NEN^  =  ^  NENf  +  NSRj  , 


(99) 


where  the  additive  noise  components  are 


MEN, 


AL 


v,det  ~ 


NEN, 


—  from  the  detector. 


from  the  amplifier. 


alias  =  ^(fimage)sl^E^ldet  +  ^EN^^^p  from  aliasing. 


(71) 

(74) 

(77) 
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(86) 


1  {'CoxT^s^l^{vyOPD)Bv)  [N 


2*"-^  '^oxVs^H^Vopd)  V  6 


—  from  quantization,  and 


NEN^^j^P  =  2n5\vB{v)\  „  from  sampling  errors. 


The  multiplicative  noise  components  are 


''v.WSL  -■ 


v(T  ypp  1  dM 


from  velocity  variations, 


and  from  optical  jitter  within  the  interferometer. 


NSRjjtt  =^  =  ^{^R?  <ya  ^^{4 }  +  <^a  ■  (98) 

Refer  to  the  preceding  subsections  for  a  full  description  of  the  variables  occurring  in  these  expressions. 
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5.  COMPARISON  TO  FILTER  WHEEL  RADIOMETER 


It  may  be  instructive  to  compare  the  noise-equivalent  spectral  radiance  of  the  FTS  to  that  of  a  filter 
wheel  radiometer  having  the  same  overall  system  parameters  [12],  For  a  filter  wheel  radiometer,  the 
signal  is 

S,=A„QT„^8v9i,B,  ,  (100) 

while  the  noise,  assumed  to  come  from  the  detector,  is 


(7^  =  NEPaet%  = 


D* 


(101) 


The  noise-equivalent  spectral  radiance  is  thus 


NEN^  fwr  - 


4^ 


(102) 


The  measurement  bandwidth,  Af,  is  inversely  proportional  to  the  integration  time  for  each  measure¬ 
ment 


Af  =  — 

iti 


Assuming  equal  dwell  times  on  each  wavenumber,  the  integration  time  is  given  by 


(103) 


ti<- 


dv 


V2-V1 


(104) 


where  an  inequality  is  shown,  because  some  of  the  total  measurement  time,  T,  should  be  spent  in 
sampling  the  background  level.  (All  variables  used  here  have  the  same  definitions  as  for  the  interferom¬ 
eter.) 

The  noise-equivalent  spectral  radiance  for  the  filter  wheel  radiometer  reduces  to 


NENy  FWR  - 


2dvT 


(105) 


For  the  Fourier  Transform  Spectrometer,  we  also  assume  that  detector  noise  is  the  dominant  term, 
as  is  generally  the  case.  It  is  given  by 


^E^v,det,FTS 


4l4a4^  it 


(106) 
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For  the  FTS,  the  noise  bandwidth  is  equal  to  the  Nyquist  frequency: 


hence 


MEN, 


A,nT„,T],5vD*l  T 
When  we  take  the  ratio  of  the  noise  levels  for  the  two  instruments,  we  find 


NEN^  fts  ^  '^D*fwr 
NEN^,fwR  Vs^*FTS 


(107) 


(108)  -  (72) 


(109) 


Here  we  have  retained  the  D*  factors  for  the  two  separate  instruments.  If  the  D*’s  in  both  cases 
are  dominated  by  thermal  (i.e.,  internal)  noise  sources,  they  would  be  equal  and  would  cancel  each  other 
in  equation  (109).  In  that  case  we  have,  assuming  that  the  beamsplitter  efficiency  is  roughly  100%, 


^£K,Frs 

NEN^,fwr 


•Det. Limit 


(110) 


where  is  the  number  of  spectral  samples,  assumed  to  be  the  same  for  both  instruments.  This  expression 
shows  the  multiplex  advantage  ascribed  to  Fourier  transform  spectrometers. 

At  the  other  extreme,  the  detector  D*’s  could  be  dominated  by  signal  photon  noise  in  both  cases. 
Then  we  have 


D* 


FTS 


^j{v2-Vi)B^ 


,  and  D*fwr  ' 


(Ilia,  b) 


where  B^  is  the  average  spectral  radiance  over  the  spectral  passband,  (vj,  V2).  The  ratio  of  the  system 
noises  becomes 


NEN.^FTS 

NEN^FWR 


^Photon. Limit 


(112) 


In  this  case,  the  Fourier  transform  spectrometer  has  a  multiplex  disadvantage,  because  each  spectral 
sample  is  affected  by  photon  noise  from  the  entire  passband.  In  addition,  there  is  the  loss  of  50%  of  the 
incident  energy  out  the  input  port  of  the  interferometer.  Equation  (112)  still  shows  the  FTS  to  be  superior 
for  measuring  a  few  strong  emission  lines.  Conversely,  the  FTS  would  be  particularly  poor  at  delineating 
narrow  absorption  features  in  a  spectrum. 
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APPENDIX  A 


APODIZATION  WITH  A  DETECTOR  FIELD-OF-VIEW  NOT  ON 

THE  OPTICAL  AXIS 

The  topic  of  self-apodization  in  an  interferometer  by  an  extended  field-of-view  was  introduced  in 
Section  2.  The  geometry  of  the  field  angles  is  shown  in  Figure  3,  which  is  repeated  here  for  convenience. 
The  simple  case,  in  which  the  circular  field-of-view  is  centered  on  the  optical  axis  (i.e.,  oLq  =  0),  is  treated 
in  Section  2.  Here  we  treat  the  more  general  case  of  ttg  ^  0. 


Figure  3.  Geometry  of  detector  acceptance  field-of-view  in  the  system  focal  plane.  Each  vector  represents  an  angle 
in  interferometer  space. 


The  starting  points  of  the  analysis  are  the  following  equations: 

■^(^)  =  2 5^,  [l  -h  G(ao ,  Po 

for  the  detected  signal,  where 

G{(Xq  ’  ^  J  cos(2mui:  cos  a)  sin  pdpd^  , 

and  the  solid  angle  of  acceptance  (included  in  the  definition  of  S^,)  is  given  by 

jQ  =  J  J^”sinpdpd^  =  2;r(l-cospo)  =  7rpo  • 


(14b) 


(16) 


(17) 
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The  spherical  law  of  cosines  gives 


cosa-cosaQCosp  +  sinaQsinpcosP  . 
When  this  is  applied  to  the  G  integral,  we  obtain 


(A-1) 


G(o;o,P))  =  ^J  cos.[Ccosp  +  SsmpcosP]smpdpdP  . 

(A-2) 

where  we  have  defined 

C  =  2kvx  cos  ttg  ^nd  S  =  2tvvx  sin  aQ  . 

(A-3a,  b) 

Expanding  the  argument  of  the  cosine  in  brackets  [...],  G  may  be  written 

^(«oPo)  =  ^(l2i  -<22) 

(A-4) 

Q\  =  \  f^°cos[Ccosp]cos[Ssinpcos)3]sinpdpdj8 

(A-5a) 

Q2~1  f^°sin[Ccosp]sin[5sinpcos)3]sinpdpd^ 

(A-5b) 

We  observe  that  in  equation  (A-5b),  the  integrand  is  antisymmetric  in  p,  so  Q2  =  0.  The  integration 
over  /?  for  Qj  yields 

Qi  =  coslCcos pJJq  (Ssin p)sin  pdp ,  (A-6) 

where  use  has  been  made  of  [7] 


cos(;c cos =  7iJq{x)  . 

In  order  to  evaluate  the  integral  in  equation  (A-6),  we  substitute 

sin  p  ,  cos  p  ,  r  2  •  2 

r  = - ,  dr  = - dp,  cosp  =  Jl-r  sin  pQ  , 

sinpo 

which  leads  to 


sinpo 


a  =2;rsin2po  [|cosrc^l-r^sin^po  Vo(-^^sinpo)-r=  /  .  dr 
L  J  Jl-r^sin^po 


We  introduce  the  approximations 


(A-7) 
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(A-8a) 


and 


1  I  2  •  2 

l  +  -r  sin  po  . 


■yjl-r^  sin^  Po  ^ 
and  neglect  sin^  Pq  relative  to  I  in  the  overall  factor  (A-8b).  We  also  make  use  of 

f2  =  2;r(l-cospo)  =  4;rsin^-^  , 


(A-8b) 


(A-9a) 


and 


TTsin^Po  =4;rsin^— cos^  — =  13  1  — sin^  — ]  =  f3  . 

2  2  V  2  J 


(A-9b) 

When  the  cosine  factor  in  (A-7)  is  expanded,  and  the  relationships  (A-8)  and  (A-9)  are  applied,  we  find 
(2i  «i3[A(C,S)  cos  C+5(C,  5) sin  C]  (A-10) 

or,  using  (A-4)  and  remembering  that  Q2  =  0, 

+  ,  (A-11) 

where 


B(C,5)  =  2j  7o(‘^''sinpo)sinf— Cr^sin^Polrdr  . 


(A- 12a) 


(A- 12b) 


Evaluation  of  the  two  r  integrals  in  equations  (A- 12)  yields  a  solution  in  terms  of  Lommel  functions 
(see  [13],  Sec.  8.8.1,  pp.  435-439).  The  Lommel  functions  are  defined  as 


s=0 


n+2s 


u_ 

w 


Jn+2s{^)  ’  and 


^«(“>vv)=  ^(-1) 

j=0 


n+2s 


W 
u 


Jn+2s{'^)  ■ 


(A- 13a) 


(A- 13b) 


The  integrals  for  A  and  B  are  found  to  be 
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(A- 14a) 


1  .  1 
cos  — M  Sin  — M 

A(C,S)  =  — j-^f/j(M,w)  +  — j-2— f/2(M.M')  ,  and 

— M 

i  o 


.  1  1 
sin  — M  cos  — M 

fi(C,5)  =  — - ^U2{u,w)  , 


where  u  and  w  are  defined  as 


We  also  use 


u  =  Cs‘m  pQ  ,  w  =  SsmpQ  . 


_ cosofosinpo  Po 

—  —  —  Sin  Pq  — - ~  — 

w  S  sintto  ccq 


Equations  (A- 14)  can  also  be  written  as 
A{C,S)  =  M{C,S) 
and 

B(C,5)  =  M(C,5)| 

where  we  define 


1  .  1  .  ■ 
cos — u  cos  i/r  +  sin  —  M  sin  ur 
2  2 


.1  1  .  ■ 

sin — u  cos  Y  -  cos — m  sin  i/r 


1  ^ 
=  M(C, 5)cos|  -u-Y 


=  M(C,  5)  sin^^  M  -  yr  I  , 


2  2 

M(C,  5)  cos  yr  = —  w)  ,  M(C,  5)  sin  yr=:  —  U2(M,w)  , 

u  u 


and 


tanyr 


^  U2{u,w) 
Ui{u,  w) 


We  observe  that 


M{C,S)  =  w)f  +  ([/2(u,  w)f  . 


Equation  (A-11)  now  becomes,  for  the  monochromatic  signal  factor  in  equation  (14), 

1 


G(ao,po}  =  M(C,S) 


cos^i  u-Y^  cos  C  +  sin|^-^  M  -  yr  I  sin  C 


,  or 


(A- 14b) 


(A-15a,b) 


(A-16) 


(A- 17a) 


(A- 17b) 


(A-18a,b) 


(A-19) 


(A-20) 


(A-21) 


34 


G{aQ,pQ)  =  M{C,S)cos(^C-^u  +  \i/j  .  (A-22) 

The  factor  M{C,S)  given  in  equation  (A-22)  is  now  seen  to  be  a  modulation  factor,  which  changes  the 
apparent  amplitude  of  the  interferogram.  The  phase  of  the  interferogram  is  also  modified  by  the  quantity 
m/2  -  ij/. 

While  equation  (A-22)  represents  an  analytical  solution  to  the  problem  of  an  off-axis  field-of-view, 
the  behavior  of  the  unusual  functions  it  contains  is  far  from  obvious.  We  need  to  examine  its  various 
components,  using  reasonable  value  ranges  for  the  constant  parameters.  First,  we  observe  from  equation 
(A-16)  that  u/w  ~  Pg/oto  ,  which,  for  GHIS,  is  much  less  than  one.  (Therefore,  the  solution  was  chosen 
over  the  V^.)  As  for  w,  the  argument  of  the  Bessel  functions  in  (A-13),  we  have 

w  s  5  sin  Po  =  Itcvx  sin  ttg  sin  pg  , 

and  both  Pg  and  are  usually  small.  (GHIS;  Pg  =  4.48  mrad,  ctg.^AX  ~  mrad.)  On  the  other  hand, 
vx  is  the  optical  delay  expressed  in  wavelengths,  and  is  related  to  the  relative  wavenumber  resolution 
(often  simply  called  “resolution”  by  spectroscopists)  as 


For  GHIS  again,  vx  is  less  than  1,000,  but  other  applications  could  require  higher  resolution.  The  maxi¬ 
mum  value  of  w  for  GHIS  is  0.671  (v  =  1150  cm'^  x  =  0.8  cm). 

With  this  regime  of  values  in  mind,  let  us  examine  the  Lommel  functions.  From  (A-13a),  we  have 


C/i(m,  w)  = 


(A-23) 


or 


Making  use  of  the  identity 

2«y„(w)  =  w/„_i(w)-t-w/„+i(w)  , 
we  find  (after  considerable  algebraic  manipulation) 


U2{u,w)  =  ^\^  )7l(w)| 


1-hH- 


\  (  u 


2  y  w 


\fu  u  T5(h') 
3vwy  lywy  •/j(w) 


(A-24) 


(A-25) 


(A-26) 


In  this  regime,  the  Bessel  functions  are  well  approximated  by  the  first  terms  in  their  power  series 
expansions: 
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(A-27) 


w 


(-ir 


2J  ;^^k\{n  +  k)\K2 


W 


2k 


nil  2 


Thus,  we  have 


yi(w)~3!ll 


y5(w)  1  ( w 


yi(H')  5!U 


,  etc. 


and  can  now  be  written 


Vw 


i-if" 

6l2 


120  U 


(A-28) 


and 


4ywJ 


i.ijr-f-iriifUj-li 

6K2J  2V2>'  120  2  12 


3U 


+  • 


(A-29) 


When  the  foregoing  expressions  are  inserted  in  equation  (A-20),  and  the  terms  higher  than  second 
order  are  dropped,  we  find 


M(C,5)  =  -yi(H-X 

w 


‘-Ki)  ^(3 


i+- 


\  f  w 


612 


(A-30) 


or  more  simply. 


M(C,5)  =  -yi(w)Jl-^ 
w  \  4o 


(A-31) 


Applying  the  power  series  for  the  Bessel  function  (to  second  order),  this  becomes 


M(C,5)- 


.2  V 


1-^ 
V 


.2^ 


1-i^- 

V  96, 


2  2 

^  w  u 

8  % 


(A-32) 


We  also  need  the  phase  angle,  y/,  from  the  definition  (A-19).  Inserting  the  approximate  expressions 
(A-28)  and  (A-29),  and  applying  a  binomial  expansion  for  the  denominator,  we  find 


u 


tany/  ~  — 


['--f-Tl 

r  2  21 

u 

,  w  u 

1-1-  —  +  — 

6K2)  i2V2y 

6[2)  _ 

~  4 

24  48 

(A-33) 


Equation  (A-22),  which  describes  the  fringe  amplitude,  finally  becomes  (with  tani/r~  y/),  after  inserting 
the  definitions  of  u,  w,  C,  and  S, 
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We  observe  from  equation  (A-34)  that  the  effective  wavenumber  in  the  interferogram  has  become  (using 
(A-9b)) 

47t 

The  foregoing  approximate  expressions  were  derived  in  order  to  elucidate  the  basic  behavior  of  the 
modulation,  M,  and  the  phase  factor  in  the  interferogram  (see  equations  (14b)  and  (A-22)).  In  practice, 
it  is  not  difficult  to  evaluate  the  actual  functions,  with,  for  example,  a  good  spreadsheet  program.  This 
has  been  done  in  WingZ  on  the  Macintosh,  for  the  particular  angles  of  the  aperture  radius  and  offset  of 
the  GHIS  spectrometer  (cCq  =  25.9  mrad,  Pq  =  4.48  mrad,  in  interferometer  space).  The  behavior  of  the 
modulation  factor  M(cCq,  p^)  is  plotted  in  Figure  A-1.  For  the  same  system,  the  residual  phase  error  in 
radians,  after  using  the  effective  wavenumber  given  in  (A-36),  is  shown  in  Figure  A-2. 


(A-36) 
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OPD  (waves) 


Figure  A-1.  The  modulation  factor  for  the  GHIS  spectrometer,  caused  by  the  offset  and  radius  of  the  detector 
aperture  (field  stop). 


0  200  400  600  800  1000 

OPD  (waves) 


Figure  A-2.  The  residual  phase  error,  in  radians,  for  the  GHIS  spectrometer,  caused  by  the  offset  and  radius  of 
the  detector  aperture  (field  stop). 


APPENDIX  B 


PERMISSIBLE  SURFACE  ERRORS  OF  COMPONENTS 


It  has  been  shown  that  the  effect  of  wavefront  errors  in  the  Michelson  interferometer  is  to  reduce 
the  fringe  contrast  from  its  ideal  value  of  one  by 

-AMs-=2n^v^(s^^  ,  (50) 

where 


(s'^  ^  =  the  mean-squared  value  of  the  wavefront  error,  or  OPD,  between  interferometer  arms,  not 
including  the  overall  optical  delay. 

Equation  (50)  can  be  used  as  a  tool  to  specify  allowable  limits  on  deviation  from  flatness  of  the 
components,  individually  or  as  a  group.  Thus,  we  need  to  determine  the  wavefront  errors  caused  by 
deviations  of  each  of  the  surfaces  in  the  interferometer.  Figure  B-1  shows  the  optical  setup,  with  labels 
for  each  of  the  surfaces. 


f 


Figure  B-1.  Typical  layout  of  an  unbalanced  Michelson  interferometer.  Wavefront  distortions  caused  by  the 
beamsplitting  surface  are  schematically  shown. 
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Some  simple  diagrams  are  helpful  for  understanding  the  OPD  effect  of  each  surface.  Figure  B-2 
shows  rays  incident  on  a  refractive  element,  whose  surface  is  displaced  by  a  distance  S  in  the  outward 
normal  direction.  Both  the  refracted  and  reflected  rays  are  shown.  By  constructing  perpendiculars  to  the 
rays  at  the  points  of  incidence  (i.e.,  parallel  to  the  wavefronts),  we  can  easily  deduce  the  increase  in 
optical  path  (OPD)  caused  by  displacement  S.  The  formulas  for  the  reflected  and  refracted  rays  are  shown 
in  the  figure.  The  OPD  in  the  reflected  rays  is  equally  applicable  to  the  Michelson  mirrors,  but  in  that 
case  we  have  cos0  =  1.  For  refracted  rays,  the  OPD  is  the  same  for  each  direction,  and  for  each  face. 

Figure  B-3  shows  the  case  of  rays  internally  reflected  from  the  displaced  surface.  Otherwise,  it  is 
the  same  as  Figure  B-2.  Notice  that  the  index  of  refraction  changes  the  OPD  effect. 


Figure  B-2.  Geometry  of  rays  refracted  and  externally  reflected  from  a  surface  with  a  positive  deviation  S. 
Construction  lines  are  drawn  parallel  to  wavefronts. 
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Figure  B-3.  Geometry  of  rays  refracted  and  internally  reflected  from  a  surface. 


For  the  interferometer  shown  in  Figure  B-1,  the  wavefront  OPD  between  the  two  arms  of  the 
interferometer  (except  for  the  overall  jc)  is  given  by 

5(r,  0)  =  2(53  (n  cos  e'  +  cos  6)  +  “  %C  +  (*^4  “  “  ^2 )(”  cos  9'  -  cos  0)]  ,  (B- 1 ) 

where 

6  =  the  angle  of  incidence  at  the  beamsplitter, 

O'  =  the  angle  of  refraction  within  the  beamsplitter  and  compensator, 

n  =  the  index  of  refraction  of  the  beamsplitter  and  compensator, 

ncos0'  =  Vn^-sin^0  ^  from  Snell’s  Law, 

S-  =  surface  deviation  of  surface  i,  measured  at  the  ray  position  (r,  0)  within  the  Lyot  stop.  Primes 
indicate  the  non-common  ray  path  intersections. 

^MS  ~  surface  deviations  of  the  mirror  on  the  beamsplitter  side  at  (r,  0),  and 

^MC  =  surface  deviations  of  the  mirror  on  the  compensator  side  at  (r,  0). 

In  equation  (B-1),  the  position  coordinates,  (r,  0),  all  refer  to  a  projected  position  within  the 
aperture.  We  notice  that,  for  small  angles  of  incidence,  the  OPD  is  twice  as  sensitive  to  beamsplitter 
errors  as  to  mirror  errors.  This  comes  about  because  a  positive  surface  deviation  causes  a  positive 
wavefront  change  in  one  leg,  and  a  comparable  negative  change  in  the  other.  The  optical  path  difference 
is  roughly  twice  as  large  as  either.  Also,  it  is  the  difference  between  the  errors  on  the  two  mirrors  which 
matters. 
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Often  it  is  easier  to  measure  (and  specify)  wavefront  irregularity  in  transmission,  rather  than 
measure  surface  irregularity.  The  transmitted  wavefront  errors  of  the  compensator  and  beamsplitter, 
respectively,  are  given  by 


^C7-(r,0)  =  (Sj -t-S2)(«cos0'-cos0)  ,  and 

(B-2) 

557'(r,0)  =  (^3 -1-54  )(n  cos  0'- cos  0)  . 

(B-3) 

Using  these  transmitted  wavefront  errors,  equation  (B-1)  becomes 

S[r,  0)  =  2[5;^^5  -  +  ^ST  ~  ^CT  +  2^3  cos  0]  . 

(B-4) 

Equation  (50)  can  be  used  to  predict  the  loss  of  fringe  contrast,  if  we  know  the  exact  surface 
irregularities  of  the  components  (using  equation  B-1),  and/or  the  transmitted  wavefront  errors  (using 
equation  B-4).  However,  to  specify  the  figure  requirements,  we  must  make  some  simplifying  assump¬ 
tions,  as  follows.  Assume  that  the  figure  errors  are  random  and  uncorrelated.  Their  mean-squared  values 
should  then  be  additive.  Also,  we  shall  express  the  errors  as  peak-to-valley  (P-V)  measurements,  given 
in  waves,  at  some  stated  test  wavelength,  We  assume  the  P-V  values  to  be  related  to  the  rms  values 
as 


for  the  flatness  of  surface  i,  or  the  transmitted  wavefront  error  of  component  j,  respectively.  Then  we  find, 
using  equation  (B-4), 


.w  Jr  ( 


.2/1  \2 


(B-6) 


where 

A  =  the  (minimum)  operating  wavelength  of  the  interferometer  (=  1/v). 

Assuming  that  both  mirrors  and  both  transmissive  elements  are  made  to  the  same  specifications, 
equation  (B-6)  can  be  simplified  to 


-AM§  <k‘' 


2 

[F^-h2F3^cos^6>-l-Wj] 


(B-7) 


All  of  the  specifications  in  equation  (B-7)  must  include  mounting  distortion  and  thermal  effects,  as  well 
as  the  original  fabrication  errors  of  the  components. 

For  the  case  of  the  GHIS  interferometer,  6  =  30°,  A  >  3.7  /xm,  and  A,g^,  =  0.633  jUm.  The  beamsplitter 
and  compensator  are  KBr,  with  n(3.7  fim)  =  1.536.  To  estimate  the  wavefront  error  effects,  let  us  assume 
that  the  mirrors  are  made  to  1/4  wave  from  fabrication,  with  a  further  1/8  wave  mounting/thermal 
distortion.  Thus,  ~  3/8.  For  the  transmitted  wavefronts,  assume  1/4  wave  from  fabrication,  and  1/20 


4 
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wave  mounting/thermal  distortion.  (It  is  not  easy  to  alter  the  transmitted  wavefront  by  bending  a  win¬ 
dow.)  Thus,  Wj~  3/10.  The  flatness  of  the  splitter  surface  is  most  critical,  and  we  assume  1/8  wave  from 
fabrication  plus  1/8  wave  distortion,  to  get  ~  1/4.  Inserting  these  assumptions  into  (B-7),  we  find 
-AM^  <  9.4%. 


♦ 
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APPENDIX  C 


MICHELSON  INTERFEROMETER  WITH  RETROREFLECTOR  MIRRORS 

The  plane-mirror  Michelson  interferometer  employed  for  FT  interferometry  is  sensitive  to  the 
angular  alignment  of  its  mirrors.  A  monochromatic  interferogram  takes  the  form 

2 

^£(v,jc)  =  — Ji(k)cos(2m'A:)  , 
k 

where 

k  =  AitveR  , 

£  =  the  tilt  misalignment  of  the  mirrors, 

R  =  the  radius  of  the  Lyot  stop, 

V  =  the  wavenumber,  and 
X  =  the  optical  delay. 

Tilt  misalignment  can  be  avoided  by  the  use  of  precise  retroreflecting  optics,  such  as  comer-cube 
mirrors,  or  cat’s-eye  reflectors.  A  Michelson  interferometer  with  comer-cube  retroreflectors  is  shown 
schematically  in  Figure  C-1.  In  principle,  retroreflectors  return  the  wavefront  exactly  parallel  to  its 
original  plane.  However,  the  returned  wavefront  is  generally  inverted  in  its  plane,  about  the  center  of  the 
retroreflector.  It  is  shown  in  this  Appendix  that  this  results  in  a  sensitivity  of  the  interferometer  to 
translational  alignment  of  the  retroreflectors. 


(36) 

(37) 
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INPUT  BEAM 


Figure  C-J.  Mic kelson  interferometer  with  comer-cube  retroreflectors  (schematic). 


,  C 


h 
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Figure  C-2.  Close-up  of  a  comer-cube  retroreflector  which  has  been  displaced  from  its  on-axis  position  (shown 
in  grey)  by  distance  s.  For  rays  incident  at  angle  9  from  the  optical  axis,  the  apparent  linear  shift  of  the  vertex  is 
shown  as  5. 
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C.1  ANALYSIS  OF  RETROREFLECTOR  DISPLACEMENT 


Consider  now  the  effects  of  translational  misalignment  of  the  retroreflectors.  (We  shall  assume  that 
the  reflectors  are  optically  perfect,  i.e.,  they  introduce  no  wavefront  distortions  by  themselves.)  Figure 
C-2  is  a  close-up  diagram  showing  one  of  the  retroreflectors  displaced  a  distance  s  from  the  optical  axis, 
which  is  defined  by  the  vertex  of  the  fixed  reflector  and  the  center  of  the  Lyot  stop.  The  displacement 
leaves  the  incoming  and  outgoing  wavefronts  parallel.  Also,  it  has  no  effect  on  a  beam  parallel  to  the  axis. 
However,  for  rays  incident  at  an  angle  0  to  the  axis,  as  shown  in  the  figure,  there  is  a  change  28  in  the 
optical  path,  relative  to  the  undisplaced  reflector.  The  optical  delay  in  the  interferometer  becomes 

x'  =  xcos0-l-25  =  xcos0-l-2ssin0  ,  (C-1) 

or,  generalizing  to  three  dimensions. 


where 


x'  =  x  cos  p  4-  2s  sin  p  cos  (j) 


(C-2) 


(p,  ^)  =  polar  and  azimuthal  angles  of  rays  in  the  field  of  view  (FOV),  and 
0  =  0  in  the  displacement  direction  of  retroreflector. 

Assume  that  an  extended  source  of  radiation  fills  the  FOV.  We  need  to  do  an  integration  over  that 
FOV  to  obtain  the  whole  signal.  The  detected  signal  in  the  interferometer  is 

j(x)  =  2j^  5v[l  +  iif5(v,x)]ifv  ,  (C-3) 

where 


//^  ( V,  x)  =  -^  J  cos[2m'(x  cos  p  +  2s  sin  p  cos  0)]  sin  pdpd<j) 


1  fPo 

n- 

Po  =  the  angular  radius  of  the  field  of  view,  and 
Q  =  2;r(l  -  cos  Po )  =  4;rsin^  =  the  solid  angle  of  the  FOV. 


Expanding*  the  outer  cosine  in  (C-4),  we  have 


where 


P/T  I*  Pq 

P\=\  cos(Ccosp)cos(5sinpcos0)sinprfpd0  , 
“—7t  Jo 

P2=j  sin(  Ceos  p)sin(5  sin  p  COS  0)  sin  pdpd0  , 


(C-4) 


(C-5) 


(C-6) 


(C-7) 


*Readers  of  Appendix  A  will  recognize  that  the  integral  in  equation  (C-4)  is  identical  in  form  to  that 
in  equation  (A- 2),  and  anticipate  that  the  solution  will  be  in  terms  of  Lommel  functions. 
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and  the  expressions  have  been  simplified  by  defining  C  s  2m’x  and  5  =  Attvs.  We  observe  that  since  the 
integrand  in  equation  (C-7)  is  antisymmetric  in  we  have  =  0.  When  we  make  use  of 


cos[x cos (j))d<l)  =  7!Jq{x)  , 


the  (j)  integration  in  yields 

Pj  =  2;r  cos(Ccosp)yo(>S'sinp)sin  pdp  . 

Jo 

Proceeding  as  in  Appendix  A,  make  the  substitution 

I  2  ^  2  1  2  2 

smpsrsinpo  ,  cosp  =  iJl-r  sin  pQ~l-  —  r  sin  Pq  , 
and  the  approximation 

2 

,  rdr  sin  pg  .  .  2 

sm  pdp=  .  =  rdr sin'^  Po  . 

yl-r  sin^po 


Thus  equation  (C-8)  becomes 


f 

C 

1  1  2  •  2 
l--r  sin  Po 

Jo  1 

1 

L  2  y 

yo(5rsinpo)r  rfr  . 


Defining  the  symbols 


2 

M  s  Csin  Po  and  w  s  5sin  po  , 
and  expanding  the  cosine,  we  have 

Pj  =;rsin^Po[A(M,H')cosC+fi(M,w)sinC]  , 


where 


fl  2^ 

cos 

Jo 

U  j 

■  f 

'1  2^ 

sm 

—  ur 

Jo  1 

a  J 

jQ{wr)rdr  ,  and 


Since  we  have 


;rsin^Po~i2  , 

equations  (C-10)  and  (C-5)  combine  (since  P2  =  0)  to  give 


(C-8) 


(C-9) 


(C-10) 


(C-lla) 

s. 

(C-llb) 


(C-12) 
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~  y4(M,  w)  cos  C  +  B(m,w)  sin  C  .  (C-13) 

Once  again,  the  two  r  integrals  in  equations  (C-11)  yield  a  solution  in  terms  of  Lommel  functions 
[13]  (Section  8.8.1.,  pp.  435-439).  The  Lommel  functions  are  defined  as 


°°  /  \n-l-2i 

Jn+lsi^)  .  and 

.=0 


(C-14a) 


5=0 


\«+25 


(“>  w)  =  S  “]  -^n+ls  ■ 


(C-14b) 


So  far,  the  solution  has  paralleled  that  in  Appendix  A;  but  here  we  are  interested  in  the  case  where 
the  retroreflector  displacement,  s,  is  small,  but  the  OPD  may  be  large,  such  that  w/u  <  1.  Thus  we  need 
a  solution  in  terms  of  the  series.  From  Bom  and  Wolf  [13],  we  find 


2  sm— w 
2  .  w  2 


2m 

1 

—  M 

2 

1 

cos  — M 
2 

2m 

1 

—  M 

2 

1 

cos  — M 
2 


A(m,w)  =  — sin - 1 — j^^Vq{u,w) - ,  ■  {u, w)  ,  and 

jj  2ij  t  1 


-u 


.  1 
sm  — M 
2 


(C-15a) 


B{u,  w)  =  —  cos - Va  (m,  w) - Vi  (m,  w) 

u  2u  ^  ^ 


—u 

2 


(C-15b) 


Inserting  these  solutions  in  (C-13)  and  collecting  trigonometric  factors,  we  find 


Hs{v,x)  =  - 
u 


sin 


V 


-  sin|^C  -  w)  -  cos|^C  -  j  Vj  (m,  w) 


(C-16) 


Though  it  is  a  solution  in  closed  form,  the  result  in  (C-16)  needs  further  reduction  before  it  can 
easily  be  compared  with  the  more  usual  expressions.  We  may  define  intermediate  symbols 

2 

u  u  w 

Y  =  C - ,  and  Z  =  -  + -  , 

2  2  2m 

and  (C-16)  becomes 


(v,  x) «  —  [sin(y  +  Z)-  sin  YVq  (m,  w)  -  cos  YVi  (m,  w)]  . 


(C-17) 


Next  we  have 


2 

Hg{v,x)~  —  |sin  F[cos  Z  -  Vq  (m,  w)]  +  cos  y[sin  Z  -  Vj  (m,  w)]}  , 


(C-18) 
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which  can  be  written 


Hg (v, x)  ~  M(u,  w)[sin  Fsin  yr  +  cos  Y cos  yr]  =  M(u,  w)cos(K  -  y/)  , 
where  we  define  (omitting  the  explicit  Lommel  function  arguments  for  brevity), 

M{u,  w)  =  —  Ji  +  Vq  +V^  -  2Vq  cos  Z  -  2Vj  sin  Z  ,  and 
u 


(C-19) 


(C-20a) 


cos  Z-Vo 

■  <C-20b) 

The  foregoing  solution  is  equivalent  to  that  given  by  Murty  [14]  for  the  same  case  of  w  <  m. 

At  this  point,  we  return  to  the  definitions  of  u,  w,  Y,  Z,  etc.: 


•  2  2 
M  =  Csin  pQ=27ivxs\n  Pq  , 

wsSsin  po=4m'5sinpo  , 

K  =  — sin^  Pol  =  2;rvjcfl 1  ,  and 


(C-21a) 

(C-21b) 


(C-21c) 


Z  =  27Cvx  — sin^Po+2f— 1  =2itvx  —  +  2(— 


(C-21d) 


In  these  expressions,  we  have  made  use  of 

.  2  ^ 

Sin  po=—  . 

71 

Equation  (C-19)  is  seen  to  be 


We  recognize  that  (C-22)  is  very  similar  in  form  to  equation  (24)  in  Section  3.2,  namely  a  fringe  cosine, 
modulated  by  a  slowly-varying  factor,  M.  The  argument  of  the  cosine  is  modified  with  the  same  factor 
as  in  (24),  but  an  additional  phase  shift  (y/)  appears. 

The  modulation  factor  M  and  phase  shift  yr  have  been  evaluated  for  the  case  of  Pq  =  1°.  and 
vx  =  1,000  waves.  The  results  are  plotted  as  a  function  of  vs  (the  displacement  in  waves)  in  Figures 
C-3  and  C-4. 


w 
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01  23456789  10 

Retroreflector  Displacement  (waves) 

Figure  C-3.  Modulation  factor  M(u,  w),  vs.  retroreflector  displacement,  vs  (in  waves),  for  the  case  ofFOV  radius 
Pq  =  1°,  and  the  interferometer  OPD  vx  =  1,000  waves  (u  =  1.91). 


Figure  C-4.  Phase  shift  yfu,  w)  vs.  retroreflector  displacement,  vs  (in  waves),  for  the  same  case  as  in  Figure 
C-3. 
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As  a  check  on  the  foregoing  result,  we  observe  that  for  j  =  0,  then  w  =  0,  Vq(u,0)  =  1,  and 
Vj(m,0)  =  0.  Equation  (C-16)  becomes 


//^  ( V,  ;c)  ~  —  [sin(  y  +  Z)  -  sin  y  ]  =  —  cos[  y  +  —  zlsinf— Z] 
u  u  \  2  J  y2  J 


or 


(C-23) 

(C-24) 


Inserting  the  definitions  above, 
H^{v,  x)  ~  sine 


f 

1 

f 

r 

vx 

—  cos 

27WX 

1 - 

V  27cJ  ' 

K 

4;rjy 

(C-25) 

which  is  identical  to  equation  (24),  and  represents  the  apodization  caused  by  the  extended  field  of  view. 

C.2  APPROXIMATE  EXPRESSIONS  FOR  THE  FRINGE  MODULATION 

To  clarify  the  behavior  of  the  amplitude  modulation  and  phase  factors,  we  may  express  them  as 
power  series  expansions,  valid  for  small  s.  We  use  the  series  for  the  Bessel  functions 

*  1  (wY 


Thus  to  first  order  in  (w/2)^, 


‘:ik\{n  +  kyX2)  n\K2 


k=0 


=  ,  and 

Inserting  these  approximations  in  (C-20a),  we  find 


(C-26a,b) 


M{u,w)-=- 

u  ^ 


2[1  -  cos  Z] 


'-'1 


-If^l  sinZ 

M  V  2 , 


(C-27) 


Applying  the  identities 

.  2  2  1  “  cos  Z 


sm  —  =  ■ 

2  2 


,  Z  sin  Z 
and  cot— =  - 


2  1  -  cos  Z 


(C-28) 


we  find 


M{u,  w)  =  —  sin  — 
u  2 


,  , wYL  2  Z 

1  -  i  —  1  +  —  cot  — 

2  y  V  u  2 


1 

2  4  .  Z 
=:  — sm  — 

u  2 


\  f  w 


1  —  —  1  +  —  cot  — 

2  V  2  y  V  M  2 


(C-29) 
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From  the  definition  of  Z, 


.  Z  .  u  u  . 

sin  — =  sin  — cos - l-cos— sm  — 

2  4  4m  4  4m 

or,  since  w  is  small  (though  u  may  be  appreciable), 


(C-30) 


.  z 

.  M 

2 

w 

M 

.  M 

1  M 

\2 

1  M 

sm —  = 

=  Sin  —  + 

COS  — 

=  sm  — 

1-I-- 

— 

cot  — 

2 

4 

4m 

4 

4 

M  ' 

.2, 

;  4 

(C-31) 


Making  the  substitution  in  (C-29),  we  get 


M(m,  w)  ~  —  sin  — 
M  4 


’iA( 

^  w' 

s2  1 

)  M 

cot— 

1- 

--( 

-Tf 

M 

j  4 

2i 

^  2  Z 
1  +  —  cot  — 
I  M  2 


(C-32) 


or 


M(m,  w)  »  —  sin  — 
M  4 


1  f  1  f  u  Z 

1 —  +—  —  cot cot  — 

2K2J  u{2)  14  2 


(C-33) 


From  the  identity 


we  find  that 


cot  A  -  cot  fi  = 


_  sin(fi- A) 


sin  A  sin  fi 


sm 


M  Z 

cot  —  cot  —  = 
4  2 


1  f  w 


M  \  2 


.  M  . 

sm  —  sml 
4 


.2^ 


u  w 

-4- - 

y^4  4m  y 


.  2  M 
sm  — 
4 


(C-34) 


Hence,  the  third  term  in  (C-33)  is  second  order  in  (w/2)^,  and  we  obtain 


For  the  phase  shift,  y/,  we  have  from  equation  (C-20b),  to  first  order  in  (>v/2)^. 


tani/r  = 


cosZ-l  +  f- 

U 


sinZ- 


2fw 


M  V  2 


(C-36) 
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which  may  be  expressed  as 


or 


or  also 


tan  yr  ~  ■ 


cos  Z  - 1 


sinZ 


1  +  1-1  - ^ - 

2  )  cos  Z  —  1 


tan  yr  =  - 


1 


2  1 
-  +  — 


cos  Z  - 1  M  sin  Z 


u\2  )  sin  Z 


Z 

tan  — 
2 


n-l 


(C-37) 


(C-38) 


tanyr : 

Let  us  assume  that 


*1+2 

1  ( 

u 

sinzl 

a) 

Z  1  (w 

tan — I - — 

2  sin  Z  V  2 


(C-39) 


tanyr  =  tan(^;i^-y^ 


tan  tan 


,  z  • 

1  +  tan  tan  Y 


To  find  we  compare  (C-39)  and  (C-40)  to  find  that 

If'' 

tan;t:^ 


w 

2y 


1  ( w 


.  ^  2fwS^  smZ'^2 
smZ  +  —  — 
m12 


Hence,  to  first  order  we  have 


1  f  w 


^  sinzU 


and 


(C-40) 


(C-41) 


(C-42) 


Z  1  fw']^ 

yr==  —  + -  — 

2  sm  Z  2  y 


_ 1 

4  VsinZ  u 


(C-43) 


Gathering  together  equations  (C-19),  (C-35),  and  (C-43),  we  have  the  approximate  expression  for 


for  small  values  of  s\ 
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r  21 

( 

/  \2') 

TJ  ^  \  4  .  M 

Hg{v,x)  =  —  sin  — 
u  4 

1-^ 

8 

cos 

C  "  1 

f  w  )  1 

(C-44) 

1  4  1 

i^sinZ  u) 

I2J  J 

We  observe  that  equation  (C-44)  is  of  almost  exactly  the  same  form  as  (C-24)  or  (24)  of  Section  3.2,  with 
the  addition  of  an  amplitude  modulation  factor  and  an  additional  phase  shift.  The  factor 

2 

l--y  =  l-2;r^(v5f  sin^Pb  (C-45) 

represents  the  loss  of  fringe  contrast  caused  by  the  retroreflector  displacement  alone. 

C.3  RETROREFLECTORS  VS  PLANE  MIRRORS 


An  obvious  question  arises  as  to  whether  better  interferometer  performance  can  be  obtained  with 
plane  mirrors  or  with  retroreflectors,  given  the  level  of  alignment  precision  that  is  attainable.  For  the 
plane-mirror  interferometer,  the  expression  for  loss  of  fringe  contrast  corresponding  to  (C-45)  is 


1  fk 


-JAk)  =  \--  -  -t 
k  ’  2\2  ^ 


'12I2 


1  -  2{7weR)^  , 


(38) 


where  R  is  the  radius  of  the  Lyot  stop,  and  e  is  the  mirror  tilt  error.  If  we  define  the  respective  losses 
of  fringe  contrast  as 


4efro=2;r2M^sin^/q),  and  Api^„^=2{7tveRf 


and  take  their  ratio,  we  find 


\etro  _r  ^sinpo  Y 
^plane  ^  ^  ' 


(C-46) 


For  any  interferometer  with  given  R  and  Pq,  equation  (C-46)  expresses  the  relative  sensitivity  to  loss  of 
fringe  contrast  caused  by  displacement  error  for  retroreflectors  vs  tilt  error  for  flat  mirrors. 


To  put  the  loss  ratio  in  perspective,  let  us  assume  that  parts  can  be  positioned  with  a  linear  precision 
of  s,  and  with  an  angular  precision  of  s/D,  where  D  is  a  characteristic  size.  Let  us  take  D  to  be  the 
diameter  of  the  Lyot  stop.  In  this  case,  the  plane  mirrors  contrast  loss  is 


^plane=^(^^f  ’ 
and  the  ratio  of  losses  becomes 

=  4  sin  Vo 

^plane 


(C-47) 


(C-48) 
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We  can  conclude  that  <  ^piane  Po  words,  the  retroreflector  misalignment  will 

cause  less  loss  of  fringe  contrast  than  plane-mirror  misalignment  under  these  assumptions,  as  long  as  the 
field-of-view  diameter  is  less  than  60°.  Since  this  is  almost  certain  to  be  true,  we  are  led  to  believe  that 
it  is  easier  to  build  a  high-performance  interferometer  with  retroreflectors,  rather  than  plane  mirrors. 

The  foregoing  conclusion  must  be  tempered  with  the  fact  that  it  is  much  more  difficult  to  fabricate 
an  ideal  retroreflector  than  an  ideal  plane  mirror.  A  complete  trade  study  must  also  take  into  account  the 
wavefront  errors  likely  to  be  introduced  by  the  different  components.  Furthermore,  the  ratio  of  the  linear 
and  angular  positioning  errors  may  be  quite  different  than  assumed  above,  when  the  details  of  the 
Michel  son  mirror  scan  motion  are  accounted  for. 


t 
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